Semiclassical statistical mechanics' tools for deformed algebras 
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In order to enlarge the present arsenal of semiclassical toools we explicitly obtain here the Husirni 
distributions and Wehrl entropy within the context of deformed algebras built up on the basis of 
a new family of g— deformed coherent states, those of Quesne [J. Phys. A 35, 9213 (2002)]. We 
introduce also a generalization of the Wehrl entropy constructed with escort distributions. The two 
generalizations are investigated with emphasis on i) their behavior as a function of temperature and 
ii) the results obtained when the deformation-parameter tends to unity. 
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I. INTRODUCTION 

The semiclassical approach has had a long and distin- 
guished history and is a very important weapon in the 
physics' armory. Indeed, semiclassical approximations 
to quantum mechanics remain an indispensable tool in 
many areas of physics and chemistry. Despite the ex- 
traordinary evolution of computer technology in the last 
years, exact numerical solution of the Schrodinger equa- 
tion is still quite difficult for problems with more than a 
few degrees of freedom. Another great advantage of the 
semiclassical approximation lies in that it facilitates an 
intuitive understanding of the underlying physics, which 
is usually hidden in blind numerical solutions of the 
Schrodinger equation. Although semiclassical mechan- 
ics is as old as the quantum theory itself, the field is 
continuously evolving. There still exist many open prob- 
lems in the mathematical aspects of the approximation 
as well as in the quest for new effective ways to apply 
the approximation to various physical systems (see, for 
instance, [J, and references therein). 

In a different vein, applications of the so-called g-calculus 
to statistical mechanics have accrued increasing inter- 
est lately Q. This g-calculus 4] has its origin in the 
g-deformed harmonic oscillator theory, which, in turn, 
is based on the construction of a SUq{2) algebra of 
g-deformed commutation or anti-commutation relations 
between creation and annihilation operators [E)S0| ■ The 
above mentioned applications also employ "deformed in- 
formation measures" (DIM) that have been applied to 
different scientific disciplines (see, for example, H, H, Q 
and references therein). DIMs were introduced long ago 
in the cybernetic-information communities by Harvda- 
Charvat [l^l and Vadja [Til in 1967-68, being rediscov- 
ered by Daroczy in 1970 [12l | with several echoes mostly 
in the field of image processing. For a historic sum- 
mary and the pertinent references see Ref. (isj . In 
astronomy, physics, economics, biology, etc., these de- 
formed information measures are often called <?-entropies 



since 1988 



In this work we are concerned with semiclassical statistical 
physics' problems and wish to add tools to the semiclas- 
sical armory by, in particular, investigating putative de- 
formed extensions of two of its most important quantities, 
namely, Husimi distributions (HD) and Wehrl entropies 
(WE) . We will work within the context of deformed alge- 
bras built up with the new family of g-deformed coherent 
states introduced recently by Quesne [141] . analyzing the 
main HD and WE properties. 



Our attention will be focused on the thermal description 
of the harmonic oscillator (HO) (together with its phase- 
space delocalization as temperature grows) , in the under- 
standing that the HO is, of course, much more than a mere 
example, since in addition to the extensive use Glauber 
states in molecular physics and chemistry (isl . [l6| , nowa- 
days the HO is of particular interest for the dynamics of 
bosonic or fermionic atoms contained in magnetic traps 
[13, dl, as well as for any system that exhibits an equidis- 
tant level spacing in the vicinity of the ground state, like 
nuclei or Luttinger liquids. 



The paper is organized as follows: in Section [IT] we in- 
troduce the behavior of the HD and the WE for the HO, 
studying the fluctuations in thermal equilibrium while, 
in Section IIIIl the new g-deformed coherent states are 
introduced. We study the behavior of our generaliza- 
tions of both the HD and the WE in Sections IIVI and 
rvl respectively. Generalized WE's bounds are analyzed 
in Section IVll while, in Section rVIIi we deal with escort 
distributions, which allow us to build up an alternative 
generalization of the WE. Finally, some conclusions are 
drawn in Section [VIIII 
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II. SEMICLASSICAL DISTRIBUTION IN 
PHASE-SPACE 



0, 1, . . . In this situation the analytic expression for the 
HD and the WE were obtained in Ref. UM 



Wehrl's entropy W is a very useful measure of localization 
in p hase-s p ace [191] . It is built up using coherent states 
[Us El and constitutes a powerful tool in statistical 
physics. The pertinent definition reads 



Kz) = (1 



W ^l- ln(l - e-'^'"^) 



(6) 



(7) 



W 



dx dp 

2TTh 



fi{x,p) lnn{x,p), 



(1) 



where fi{x,p) = (^IpI^) is a "semi-classical" phase-space 
distribution function associated to the density matrix p 
of the system 15, 20]. Coherent states are eigenstates of 
the annihilation operator a, i.e., satisfies a\z) — z\z). 
The distribution p{x,p) is normalized in the fashion 



(dx dp/2Trh) ijl{x,p) = 1, 



(2) 



and it is often referred to as the Husimi distribution [22] ■ 
The last two equations clearly indicate that the Wehrl 
entropy is simply the "classical entropy" (|T]) of a Wigner- 
distribution. Indeed, is a Wigner-distribution 

Dw smeared over an h sized region of phase space [2ll |. 
The smearing renders iJ,{x,p) a positive function, even if 
Dw does not have such a character. The semi-classical 
Husimi probability distribution refers to a special type of 
probability: that for simultaneous but approximate loca- 
tion of position and momentum in phase space [2l| . The 
uncertainty principle manifests itself through the inequal- 
ity 1 < W, which was first conjectured by Wehrl IQ] and 
later proved by Lieb [l^] . 

The usual treatment of equilibrium in statistical mechan- 
ics makes use of the Gibbs's canonical distribution, whose 
associated, "thermal" density matrix is given by 



P 



z- 



(3) 



with Z = TT{e-^") the partition function, /3 = l/fcsT 
the inverse temperature T, and ks the Boltzmann con- 
stant. In order to conveniently write down an expres- 
sion for W consider an arbitrary Hamiltonian H of eigen- 
energies En and eigenstates \n) [n stands for a collection 
of all the pertinent quantum numbers required to label 
the states). One can always write pl| 



H{x,p) 



-0E„ 



(4) 



A useful route to W starts then with Eq. ^ and con- 
tinues with Eq. ([T]). In the special case of the harmonic 
oscillator the coherent states are of the form [201 



When T ^ 0, the entropy takes its minimum value 
W = 1, expressing purely quantum fluctuations. On the 
other hand when T —^ oo, the entropy tends to the value 
— hi{phuj) which expresses purely thermal fluctuations. 



III. QUESNE'S NEW «j-DEFORMED 
COHERENT STATES 

Quesne advanced in Ref. [3] a new family of harmonic 
oscillator physical states, labelled by < g < 1 and 

z € C, 



E 



vW? 



(8) 



where \n) = (n!) ^/^(a^)"|0) is the n-boson state and 
one introduces quasi-factorials [13] 



[n]g[n-l]g...[l]g 
1 



if n = 1,2,3. 
if n = 



with 



1 - (7"" 



(9) 



{n}q is called the "g-basic number" and generates its 
own factorial (the g-factorial) {n}q\ [13]. The g-factorial 
can also be written in terms of the g-gamma func- 
tion r,(x) [13 



]n],! =g^^^r,(n + l) =g' 



"M,!, (10) 



that, in the limit g ^ 1, yields [n]q, [n]q\, and Tq{n) tend- 
ing to n, n\, and T{n), respectively. The states \z)q allow 
one to build up g-deformed coherent states. They will 
be acceptable generalized coherent states if three basic 
mathematical properties are verified [IB]. The states \z)q 
should be: i) normalizable, ii) continuous in the z-label. 
Additionally, iii) one must ascertain the existence of a 
resolution of unity with a positive weight function. The 
normalization condition, ^{z\z)q = 1 leads to 



71=0 



(5) 



where \n) are a complete orthonormal set of eigenstates 
and whose spectrum of energy is En = {n + l/2)fiw, n = 



N,i\z\')^Y.^^Eq[{l^q)q\z\^^ 



(11) 



where Eq{x) = nfe°=o(l"'"'?'^-^) of the so-called Jack- 

son's g-exponentials introduced in 1909 ^26| such that 
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livaq^i Eq[{l — q)x] = e^. Since A^q(|zp) is equal to 
— (a well defined function in < g < 1), the 
new states are normalizable on the whole complex plane. 
On the other hand, the states \z)q are always continuous 
in z. Finally, for the resolution of unity one needs 



d2z|z),i^,(|z|2),(z|^^|n)(n|=/, (12) 

n=0 

with a weight function i^q(|zp) that can be obtained us- 
ing the expressions for \z)q and its conjugate and then 
performing the integral by recourse to the g-analogue of 
the Euler gamma integral Q • Quesne finds 



1 ■ t=i 

0.8 



/J, ( I J I ) 
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FIG. 1: tJ.q{z) vs. \z\ for two different normalized temperatures 
t = T/hw, and several values of q. 




l-q Eq[{l - q)q\z\^] 
TTlng-i Eq[{l - q)\z\^] ' 



(13) 



It is seen that, in the limit q ^ 1^, we have X^dzp) 
K{\z\'^) = I/tt, corresponding to the weights for the con- 
ventional coherent states of the harmonic oscillator. This 
entails that we indeed have at hand new g-deformed co- 
herent states that fulfill the standard properties. In the 
next section we start presenting the results of this com- 
munication, using these Quesne-coherent states to gener- 
alize our g-Husimi distribution. 



We see in Fig. [T]that a g-HD's height depends only upon 
the temperature (not on g!). g-Deformation affects only 
shapes. Thus, 



< ^J-q{z) < 1, 



(18) 



and ^q{z) remains a legitimate semi-classical probability 
distribution (unlike Wigner's one). 



DEFORMED WEHRL ENTROPIES 



IV. THE g-HUSIMI DISTRIBUTION TO < g < 1 

We define the g-Husimi distribution (g-HD) in the rather 
"natural" fashion 



Hq{z) = q{z\p\z)q, 



(14) 



using Quesne's HO-g-coherent states ([5]). From these it is 
easy to find an analytic expression for our new g-Husimi 
distributions 



,,(.). (1-e^^-) Ml 



Eq[{l-q)q\z\^ 
It is important to remark that 



lim^g(z) = ^l{z) 

q^l 



(15) 



(16) 



where fi{z) is given by We have numerically ascer- 
tained that our g-Husimi distribution is normalized using 
the same weight function iiTgdzp) that Quesne employed 
for the resolution of unity of her g-coherent states. 



d^zKq{\z\'')^lq{z) = 1. 



(17) 



In order to introduce deformed Wehrl entropies we start 
by reminding the reader of the concept of g-deformed 
measures, that read [291 



Sq = ~kB / dnfixY Inqfix), 



(19) 



where ks is the constant of Boltzmann's, /(x) a normal- 
ized probability density, and one defines the g-logarithm 
function in the fashion [2^ 

In^ .T = (1 - x^-'^)/{q - 1) with x> 0; g G 3?, (20) 

where g is called the deformation index (the signature 
of the "deformed" nature of the measure). For q = 1 
we reobtain the ordinary logarithm and the logarithmic 
Shannon measure. The measure (|19p can also be gener- 
ated using the g— calculus introduced by Jackson [26| in 
the following form 



Sq = ~kBDi^^g{a) 



a=l' 
(<?) 



(21) 



where g{a) — J dxf{x)°' and Dx is Jackson's g- 
derivative operator 



D 



(q) ^ /t(gx) - h{x) 
qx — X 



(22) 



which reduces to the ordinary Leibnitz-one d/dx when 
the parameter g goes to unity ^27i] . 
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Accordingly, the natural definition of a deformed Wehrl 
entropy reads as follows 



Wr, 



(fzKq{\z\^)llq{z) In Aiq(z), 



(23) 



and, inserting the explicit HO-form (llSp into the above 
expression, we find 



Wq = /(g) - ln(l ~ e-^^) = f{q) - 1 + W^, (24) 



where 



/(?) = - / Kq(PzHq\TL 



(25) 



It easy to check that in the limit q ^ 1 one has f{q) 1 
and Eq. leads to the usual Wehrl entropy. The q- 

Wehrl entropy can not be obtained in analytic fashion 
and needs numerical evaluation. In Fig. [2] we plot Wq 
vs. t [t — T jhio) for several values of the parameter q. 
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the difference between Tsallis and Shannon entropies be- 
comes weaker than either at the quantal or the classical 
levels. 



VI. WEHRL ENTROPY BOUNDS 

Using the definition of the q-WE given by Eq. ((24| we can 
easily investigate possible bounds for the Wehrl entropy, 
i.e., for our localization power. We see that, when T — ^ 0, 
one has 



where 

9(9) 
while, for T 



d\. 



\^{Eq{{\-q)q\z\^\) 
' Eq\{\-q)\zY\ ' 



00, 



— ln(/3fia;) = cx). 



(26) 



(27) 



(28) 



having a lower bound only. Since, obviously, Wq is a 
monotonously growing function of T we can state that 



Aq < Wq, 



(29) 



which constitutes the new deformed Lieb-relation. It is 
easy to check the limit q ^ 1: 



lim Ao 



1, 



(30) 



and we reobtain the known Licb bound of the standard 
Wehrl entropy 1 <W. In Tabic U we illustrate the behav- 
ior of the function Aq vs. q. We note that the function 
reproduces the values of the Wehrl entropy at T = 0, 
as expected. 



FIG. 2: Wgi^lg 



t (t — T/huj), for values ranging down- 



wards from q = 0.95 q — 0.50. 
{q = 1) is also represented. 



The standard standard WE 



We see from psp and comparison with the conven- 



tional HO- Wehrl entropy W 



ln(l 



) that 



deformation merely entails a constant-in-z additive func- 
tion f{q), depending solely on q, as depicted in Fig. [21 so 
that Wq behaves, as a function of deformation, exactly 
like the g-distributions of the previous Sections. At low 
enough temperatures our g-coherent states minimize the 
g-entropy, i.e., Wq ~ 1, yielding maximal localization. 
Delocalization grows with T, of course. 

At this stage we are in a position to draw an important 
consequence. In general, deformed (or Tsallis') entropies 
are quite different objects as compared with Shannon's 
one, save for q close to unity. This is not the case for 
Wehrl entropies! g-deformation becomes in this case a 
smooth shape-deformation. At the semiclassical level 



TABLE I: Aq vs. 
_g Aq_ 

0.50 1.5289 
0.55 1.4337 
0.60 1.3536 
0.65 1.2852 
0.70 1.2264 
0.75 1.1753 
0.80 1.1307 
0.85 1.0916 



VII. ESCORT GENERALIZATIONS 

In this Section we consider a new possible deformation of 
the entropy, which we shall call g-escort Wehrl entropy 
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(g-eWE) . It is build up from a deformation of the HD ad- 
vanced in Ref. [13], called the g-escort Husimi distribu- 
tion (g-eHD). In general, given a normalized probability 
density (PD) /(x), its associated q-escort PD F{x) reads 



3.5 



(31) 



If we replace here f{x) by the Husimi distribution one 
gets HI 

j,{\z\) = q{l - e-P^-) exp[-g|z|^(l - e-P'^-)], (32) 
where q £ (1,^2) [13]. 

Our goal is now to get the associated q-Wehrl's entropy. 
Replacing Eq. ([5^ into Eq. ([T|), we immediately find 



W(^M ^ q'~il~e-^'-^)'l,~Hqil-e-^'-n]N,, (33) 
where 



-lzl^exp[-9l^l'(l-e-^, (34) 



and 



7,(1^1). 



(35) 



Integrating over all phase space one finds A^^ = 1, the 
normalization condition of the (7-eHD. On the other hand 
we have = l/q^(l— e"''''")^. Thus, our q-eWE, is given 

by 



M^(,)(7,) = l-ln[g(l- 



)]=I^-ln(g), (36) 



a rather interesting result that entails, for our alternative 
escort q-deformation, that it simply adds a new term to 
the entropy, given by \n{q) (vanishing, as it should, for 
9 = 1)-. 

Analyzing the behavior of the g-eWE with temperature 
we can look for its bounds and easily ascertain that 



T 



^ 
00 



1 - ln(<7) 

— hi{[3hLj) = 00, 



i.e., only a lower bound exists. In table HIl we illustrate 
how the lower bound changes with q. Fig. ^depicts W^(g) 
vs. t for several values of q. The g-cscort entropies is able 
to "perforate" the Lieb lower bound for q> 1, without 
violating the uncertainty principle, since q G (0, 



VIII. CONCLUSIONS 




FIG. 3: q-EWe vs. t (t = T/huj), for q = l.l to q = 1.4. The 
standard WE (q = I) is also represented 



TABLE II: 1 — In(g) for different values of i 

q 1 - ln(g) 

1.1 0.904690 

1.2 0.817678 

1.3 0.737636 

1.4 0.663528 



• we have advanced a ^-generalization of the Husimi 
distribution ^iq{z)^ which arises from the Quesne- 
family of g-coherent states and found that the q- 
deformation does not change the HD's property of 
being legitimate probability distributions, i.e., < 

< 1- 

• the above leads to a concomitant generalization of 
the WE Wq{iiq), whose lower bound coincides with 
the well-known Lieb one. These semiclassical q- 
entropies approach the standard one when q tends 
to unity. 

• a different generalization of the Wehrl entropy, 
IF(q)(7g), called the escort-one has also been ad- 
vanced in Section fVIII starting from the g-eHD of 
Ref. [13]. We saw that this generalization closely 
resembles the standard WE, but its lower bound 
improves on the Lieb one, allowing (at least in 
principle) for a better localization in phase-space. 
This point should be further considered in future 
research. 



Semiclassical Husimi distributions and their associated 
Wehrl entropies have been here investigated within the 
frame of deformed algebras. As a summary of our results 
we can state that: 
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